Method and device for acquiring precoding matrix

ABSTRACT

A method and device for acquiring a precoding matrix are provided by the present invention. The method for acquiring a precoding matrix includes: for each receiving end, selecting the interference vectors as the interference space basis vectors from whole interference vectors, representing the remaining interference vectors as the linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors ( 101 ); determining the currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition ( 102 ). The technical solution of the present invention can eliminate interference and do not generate Bit Error Rate (BER) floor.

CROSS-REFERENCE TO RELATED APPLICATION

This application is a continuation of International Patent Application No. PCT/CN2011/073733, filed on May 6, 2011, which claims priority to Chinese Patent Application No. 201010168681.4, filed on May 7, 2010, both of which are hereby incorporated by reference in their entireties.

FIELD OF TECHNOLOGY

The present invention relates to the field of communication technology, and particularly, to a method for acquiring a precoding matrix and a device therefore.

BACKGROUND OF THE INVENTION

The Multi Input Multi Output (MIMO) technology is a significant breakthrough of the smart antenna technology in the field of wireless mobile communication. The MIMO technology can improve the communication system capacity and the spectrum efficiency exponentially without increasing the bandwidth, and it is the key technology applied to the new generation of mobile communication system. Both the transmitting end and the receiving end of the MIMO system are provided with a plurality of antennas. During the signal transmission, the transmitting end may transmit a plurality of independent signals through different transmitting antennas in a Spatial Multiplexing (SM) mode, while the receiving end may receive signals from different receiving antennas and acquire signals transmitted by the transmitting end from the received signals.

In the prior art, the transmitting ends only share their respective channel information. Then each transmitting end centralizes interferences from all other transmitting ends to the signal spaces that can be processed by corresponding receiving ends in an Interference Alignment (IA) method based on a principle of maximizing the Signal-to-Interference and Noise Ratio (SINR) at the receiving end and a principle of maximizing the Signal-to-Leakage and Noise Ratio (SLNR) at the transmitting end.

The prior art has the following deficiency:

In the prior art, each transmitting end centralizes the interferences from all other transmitting ends to the signal spaces that can be processed by corresponding receiving ends based on the principles of SINR and SLNR maximization, and can not completely eliminating interferences. As a result, when the Signal to Noise Ratio (SNR) is lowered to a certain degree such as 40 dB or 50 dB, the Bit Error Rate (BER) is no longer reduced along with the decrease of the SNR, and a BER floor occurs.

SUMMARY OF THE INVENTION

The embodiments of the present invention provide a method for acquiring a precoding matrix and a device therefore, which are capable of eliminating signal interferences and avoiding the BER floor.

For this purpose, the embodiments of the present invention provide:

A method for acquiring a precoding matrix, comprising:

for each receiving end, selecting interference vectors as interference space basis vectors from whole interference vectors, and representing remaining interference vectors as a linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors; and

determining a currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition.

A device for acquiring a precoding matrix, comprising:

an acquisition unit configured to select, for each receiving end, interference vectors as interference space basis vectors from whole interference vectors;

a linear combination unit configured to represent remaining interference vectors as a linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors; and

a calculation unit configured to determine a currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition.

In the embodiments of the present invention, for each receiving end, the interference vectors to be aligned are all represented as a linear combination of the interference space basis vectors, thus the interferences can be completely aligned and eliminated, without a BER floor occurring.

BRIEF DESCRIPTION OF THE DRAWINGS

In order to more clearly describe the technical solutions of the embodiments of the present invention or the prior art, the drawings to be used in the descriptions of the embodiments or the prior art are briefly introduced as follows. Obviously, the following drawings just illustrate some embodiments of the present invention, and a person skilled in the art can obtain other drawings from these drawings without paying a creative effort.

FIG. 1 is a flowchart of a method for acquiring a precoding matrix provided by an embodiment of the present invention;

FIG. 2 is a schematic diagram of an Interference Alignment (IA) provided by an embodiment of the present invention;

FIG. 3 is a simulation diagram of an IA scheme provided by an embodiment of the present invention; and

FIGS. 4A and 4B are structure diagrams of a device for acquiring a precoding matrix, respectively.

DETAILED DESCRIPTION OF THE INVENTION

The technical solutions of the embodiments of the present invention will be clearly and completely described as follows with reference to the drawings. Apparently, the described embodiments are just a part of embodiments of the present invention rather than all the embodiments. Based on the embodiments of the present invention, any other embodiment obtained by a person skilled in the art without paying a creative effort will fall within the protection scope of the present invention.

Referring to FIG. 1, an embodiment of the present invention provides a method for acquiring a precoding matrix, including:

101: for each receiving end, selecting interference vectors as interference space basis vectors from whole interference vectors, and representing remaining interference vectors as a linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors.

102: determining a currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition.

Wherein, the interference vector may be represented as a product of a channel vector and a precoding vector. In that case, for the kth receiving end, it is determined in step 101 that the number of the interference space basis vectors is N−L_(K), and N−L_(K) interference vectors are selected as the interference space basis vectors from L−L_(K) interference vectors; N is the number of the receiving antennas of the kth receiving end, L_(K) is the number of the signal flows of the kth receiving end, and L is the total number of the signal flows of a plurality of receiving ends. Specifically, since the number of the receiving antennas of the kth receiving end is N and the number of the signal flows of the kth receiving end is L_(K), the kth receiving end has an L_(K)-dimensional signal space and N−L_(K)-dimensional interference space. Since the total number of the signal flows of the plurality of receiving ends is L and the number of the signal flows of the kth receiving end is L_(K), the kth receiving end totally has L−L_(K) interference vectors, from which N−L_(K) interference vectors are selected as the interference space basis vectors, thus the number of the remaining interference vectors (i.e., the interference vectors to be aligned) is (L−L_(K))−(N−L_(K))=L−N.

In that case, the process of step 102 may specifically includes:

acquiring a matrix

$\begin{bmatrix} {{\overset{\_}{H}}_{1}\alpha_{1}} \\ {{\overset{\_}{H}}_{2}\alpha_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - N})}K}\alpha_{{({L - N})}K}} \end{bmatrix} = 0$

according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquiring a first Lagrange function in an optimization method according to the matrix and the current channel condition;

wherein, α_(n)=└α_(n,0 α) _(n,1 α) _(n,2) . . . α_(n,(N-L) _(k) ₎┘, n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, α_(n,0)=−1, α_(n,1), α_(n,2), α_(n,3) and α_(n(N-LK)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors;

acquiring that Aw=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end; acquiring a second Lagrange function in an optimization method according to the fact that Aw=0 and the current channel condition; wherein, each submatrix A_(nl) in A is N×M dimensional, l≦L, L is the total number of the data flows; an A_(nl) corresponding to a remaining interference vector in the nth equation is a product of α_(n,0) and a channel matrix in the remaining interference vector, while an A_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector;

acquiring a cross-correlation function between vector a_(n) and precoding matrix w using the first and second Lagrange functions; wherein specific forms of the first and second Lagrange functions will be detailedly described in subsequent embodiments;

substituting a preset initial value of w or a preset initial value of vector a_(n) into the cross-correlation function to acquire a current value of w or a current value of vector an, then substituting the current value of w or the current value of vector a_(n) into the cross-correlation function until a preset condition is satisfied. The preset condition is that the difference between a value acquired by substituting the acquired current value of w or the acquired current value of vector a_(n) into the first or second Lagrange function this time and a value acquired by the substitution last time reaches a first preset threshold, or that the number of iterations reaches a second threshold.

Or, the interference vector may be represented as a product of a channel vector and an equilibrium vector. In that case, for the kth receiving end, it is determined that the number of the interference space basis vectors is M−L_(K), and M−L_(K) interference vectors are selected as the interference space basis vectors from L−L_(K) interference vectors; wherein, M is the number of transmitting antennas of a transmitting end pairing with the kth receiving end, L_(K) is the number of the signal flows of the kth receiving end, and L is the total number of the signal flows of a plurality of receiving ends. Specifically, since the number of the transmitting antennas of the transmitting end pairing with the kth receiving end is M and the number of the signal flows of the kth receiving end is L_(K), the transmitting end has an L_(K)-dimensional signal space and M−L_(K)-dimensional interference space. Since the total number of the signal flows of the plurality of receiving ends is L and the number of the signal flows of the kth receiving end is L_(K), the transmitting end totally has L−L_(K) interference vectors, from which M−L_(K) interference vectors are selected as the interference space basis vectors, thus the number of the remaining interference vectors (i.e., the interference vectors to be aligned) is (L−L_(K))−(M−L_(K))=L−M.

In that case, the process of step 102 may specifically includes:

acquiring a matrix

$\begin{bmatrix} {{\overset{\_}{H}}_{1}\beta_{1}} \\ {{\overset{\_}{H}}_{2}\beta_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - M})}K}\beta_{{({L - M})}K}} \end{bmatrix} = 0$

according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquiring a third Lagrange function in an optimization method according to the matrix and the current channel condition;

wherein, β_(n)=└β_(n,0) β_(n,1) β_(n,2) . . . β_(n,(M-L) _(k) ₎┘, n=(k−1)(L−M)+n′, n′=1, 2, . . . , L−N, β_(n,0)=−1, β_(n,1), β_(n,2), β_(n,3) and β_(n(M-LK)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors;

acquiring that BG=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end; and acquiring a fourth Lagrange function in an optimization method according to the fact that BG=0 and the current channel condition; wherein, each submatrix B_(nl) in B is M×N dimensional, l≦L, L is the total number of the data flows; a B_(nl) corresponding to a remaining interference vector in the nth equation is a product of β_(n,0) and a channel matrix in the remaining interference vector, while a B_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector; in which, G is an equilibrium matrix composed of equilibrium vectors;

acquiring a cross-correlation function between vector β_(n) and the equilibrium matrix is acquired using the third and fourth Lagrange functions;

substituting a preset initial value of equilibrium matrix G or a preset initial value of vector β_(n) into the cross-correlation function to acquire a current value of equilibrium matrix G or a current value of vector β_(n), then substituting the current value of equilibrium matrix G or the current value of vector β_(n) into the cross-correlation function until a preset convergence condition is satisfied. The preset convergence condition is that the difference between a value acquired by substituting the acquired current value of equilibrium matrix G or the acquired current value of vector β_(n) into the third or fourth Lagrange function this time and a value acquired by the substitution last time reaches a third preset threshold, or that the number of iterations reaches a fourth threshold.

In which, the transmitting end in the embodiments of the present invention may be a base station, and the receiving end may be a User Equipment (UE).

In the embodiments of the present invention, for each receiving end, the interference vectors to be aligned are all represented as a linear combination of the interference space basis vectors, thus the interferences can be completely aligned and eliminated, without any BER floor occurring.

I. In order to describe the idea that in step 101, the interference vectors to be aligned shall be represented as a linear combination of the interference space basis vectors, the following example is given in which each UE only receives one signal flow.

The signal received by the kth receiving end may be represented as:

$\begin{matrix} {y_{k} = {{H_{kk}W_{k}x_{k}} + {\sum\limits_{{i = 1},{i \neq k}}^{K}{H_{ki}W_{i}x_{i}}} + n_{k}}} & (1) \end{matrix}$

in which, H_(kk)W_(k)x_(k) is a received useful signal,

$\sum\limits_{{i = 1},{i \neq k}}^{K}{H_{ki}W_{i}x_{i}}$

is a received interference signal, and n_(k) is a noise vector having a variance σ². Herein H_(kk) represents a channel matrix from the kth Base Station (BS) to the kth UE, and H_(ki) represents a channel matrix from the ith BS to the kth UE, with the dimensions of N×M;

W_(k) is a precoding matrix of the kth BS, and x_(k) is a signal transmitted by the kth BS; W_(i) is a precoding matrix of the ith BS, and x_(i) is a signal transmitted by the ith BS.

In which, in Formula (1):

${W_{i}x_{i}} = {\sum\limits_{m = 1}^{L_{i}}{w_{im}x_{im}}}$

L_(i) is the number of the precoding vectors of the ith UE, namely the number of the signal flows transmitted by the ith UE; all the UEs (assuming that the number of each type of UEs is K) totally have

$L = {\sum\limits_{m = 1}^{K}L_{i}}$

signal flows; it is not hard to understand that the L precoding vectors shall be pairwise independent; w_(im) is the mth precoding vector of the ith BS, with the dimensions of M×1; and x_(im) represents a transmitted signal, which is a scalar quantity.

Firstly, the IA scheme provided by the present invention is explained in a case of L_(i)=1, i.e., each UE just has one flow, thus Formula (1) shall be converted to be:

y _(k) =H _(k1) w ₁ x ₁ +H _(k2) w ₂ x ₂ + . . . +H _(kk) w _(k) x _(k)+ . . . +n_(k)

the transmitted signal estimated according to the received signal of the kth UE is:

{circumflex over (x)} _(k) =g _(k) ^(H) y _(k) =g _(k) ^(H) H _(kk) w _(k) x _(k)+(g _(k) ^(H) H _(k1) w ₁ x ₁ +g _(k) ^(H) H _(k2) w ₂ x ₂+ . . . )+g _(k) ^(H) n _(k)  (2)

In which, g_(k) ^(H) is a 1×N-dimensional equilibrium vector for weighting. When ∥g_(k)∥²=1, the selection of g_(k) does not change the noise power, so there is an essential condition for acquiring the optimal estimation of x_(k): for any i,k=1, 2, . . . , K:

g _(k) ^(H) H _(kl) w _(i)=0, i≠k  (3)

g _(k) ^(H) H _(kk) w _(k) =c _(k) , c _(k) is a positive number;

in which, in order to ensure the maximization of the SINR, c_(k) shall be maximized to achieve the optimal detection performance.

For the kth UE, when w_(i) is determined, a_(i)=H_(ki)w_(i) may be assumed so as to satisfy g_(k) ^(H)H_(ki)w_(i)=0|_(i≠k), i=1, 2, . . . , K, thus g_(k) ^(H)a_(i)=0∥_(i≠k), wherein a_(i).=(α_(i1), α_(i2), . . . , α_(iN))^(T), and a linear equation system composed of K−1 equations can be acquired as follows:

$\begin{matrix} \left\{ \begin{matrix} {{{a_{11}g_{k\; 1}^{*}} + {a_{12}g_{k\; 2}^{*}} + \ldots + {a_{1N}g_{kN}^{*}}} = 0} \\ \ldots \\ {{{a_{{k - 1},1}g_{k\; 1}^{*}} + {a_{{k - 1},2}g_{k\; 2}^{*}} + \ldots + {a_{{k - 1},N}g_{kN}^{*}}} = 0} \\ {{{a_{{k + 1},1}g_{k\; 1}^{*}} + {a_{{k + 1},2}g_{k\; 2}^{*}} + \ldots + {a_{{k + 1},N}g_{kN}^{*}}} = 0} \\ \ldots \\ {{{a_{K\; 1}g_{k\; 1}^{*}} + {a_{K\; 2}g_{k\; 2}^{*}} + \ldots + {a_{KN}g_{kN}^{*}}} = 0} \end{matrix} \right. & (4) \end{matrix}$

In order to ensure untrivial solutions of Formula (4), it is required that among the K−1 equations, the coefficient vectors of at least L−N interference zero-forcing equations are linear combination of the coefficient vectors of N−1 interference zero-forcing equations, which is the IA constraint condition. In which, the interference zero-forcing equations refer to equations with zeros on the right of the equal signs in Formula (4). In case one UE only receives one signal flow,

$\begin{matrix} {{H_{k,i_{m}}w_{i_{m}}} = {\sum\limits_{n = 1}^{N - 1}{\alpha_{kn}H_{k,i_{n}}w_{i_{n}}}}} & (5) \end{matrix}$

In which, k=1, 2, . . . , K, i_(n)|_(n=1, 2, . . . , N-1)εS_(k1), i_(m)|_(m=1, 2, . . . , K-N)εS_(k2), α_(kn) represents a combination coefficient, S_(k1) represents a set of interference space basis vectors and the set includes N−1 elements, S_(k2) represents a set of interference vectors to be aligned and the set includes K−N elements, and S_(k1)∪S_(k2)∪{k}={1, 2, . . . , K}.

For the convenience of description, the above scheme adopts a single flow as an example, i.e., each UE just receives only one path of signals, namely each UE has one signal flow. In case the kth UE has L_(K) signal flows, the total number of the interference vectors is L−L_(K) and the number of the interference space basis vectors is N−L_(K) for the UE, thus N−L_(K) interference vectors are selected as the interference space basis vectors from the whole interference vectors, and the remaining interference vectors are represented as the linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors to be aligned and the number thereof is (L−L_(K))−(N−L_(K))=L−N. In which, N is the number of the receiving antennas of the kth UE, and L is the total number of the signal flows of all the UEs.

The following example is given to make the IA scheme provided by the embodiment of the present invention be clearer. FIG. 2 is a schematic diagram of an IA provided by an embodiment of the present invention. Assuming that M=6, N=3 and K=5, i.e., there are five BSs and five UEs, wherein each UE has three receiving antennas and just receives one signal flow. For each UE, since the number of the receiving antennas is 3, the received signal vector has three dimensions, wherein one dimension is occupied by the useful signal, and two dimensions are reserved. In order to realize the IA, any three of the four interference vectors shall be linearly correlated. For this purpose, two of the four interference vectors are arbitrarily selected as the interference space basis vectors, and the two remaining interference vectors (i.e., the interference vectors to be aligned) are represented as the linear combination of the interference space basis vectors.

Referring to FIG. 2, for the first UE as illustrated in the first row, H₁₄w₄ is an interference vector to be aligned, and H₁₄w₄ is a linear combination of H₁₂w₂ and H₁₃w₃ and may be represented as:

H ₁₄ W ₄=α₁₁ H ₁₂ w ₂+α₁₂ H ₁₃ w ₃  (6)

Similarly, H₁₅w₅ is a second interference vector to be aligned, and H₁₅w₅ is a linear combination of H₁₂w₂ and H₁₃w₃, and may be represented as:

H ₁₅ w ₅=α₁₃ H ₁₂ w ₂+α₁₄ H ₁₃ w ₃  (7)

For the first UE, the number of the aligned interference vectors is L−N=2, and the number of the remaining interference vectors is N−1=2. Similarly, for the second, third, fourth and fifth UEs, multiple equations like (6) and (7) can be acquired to form an equation system.

II. Next, the situation of each UE receiving a single signal flow is extended to a situation of each UE receiving multiple signal flows. Assuming that the number of the signal flows received by the kth UE is Lk (Lk≧1), thus the signal flows of the UE occupy an L_(k)-dimensional space (the L_(k)-dimensional space is called as signal space), and the interference vectors shall be compressed into remaining N−L_(k)-dimensional space (the L_(k)-dimensional space is called as interference space), i.e.:

H _(ki) w _(im) , m=1, 2, . . . , L _(i) , i=1, 2, . . . , K, i≠k

In other words, for the kth UE, L−L_(k) interference vectors must be aligned into the N−L_(k) dimensional interference space. Assuming that π_(k) represents a set of N−L_(k) vectors selected as the interference space basis vectors from L−L_(k) interference vectors, and π_(k)′ represents a set of the remaining (L−L_(k))−(N−L_(k))=L−N interference vectors to be aligned. For example, if K=3, M=10, N=15, L=18 and L1=L2=L3=6, for each UE, the number of the interference vectors is L−L_(k)=12, and the interference space is N−L_(k)=9 dimensional, i.e., the number of the interference space basis vectors is 9, thus π_(k) and π_(k)′ for each UE may be represented as:

π₁ ={w ₂₄ ,w ₂₅ ,w ₂₆ ,w ₃₁ , . . . , w ₃₆}, π₁ ={w ₂₁ ,w ₂₂ ,w ₂₃}

π₂ ={w ₁₄ w ₁₅ ,w ₁₆ ,w ₃₁ , . . . , w ₃₆}, π₂ ={w ₁₁ ,w ₁₂ ,w ₁₃}

π₃ ={w ₁₄ ,w ₁₅ ,w ₁₆ ,w ₂₁ , . . . , w ₂₆}, π₃ ={w ₁₁ ,w ₁₂ ,w ₁₃}

in such a representation method, the IA may be represented as:

{H _(kj) w _(jm′) :w _(jm′)επ′_(k)}⊂span[H _(ki) w _(im) :w _(im)επ_(k) ], k=1,2, . . . ,K  (8)

Formula (8) indicates that the vectors in π_(k)′ are represented as a linear combination of the vectors in π_(k), thus Formula (8) is equivalent to the following vector equation:

$\begin{matrix} {{{H_{kj}w_{{jm}^{\prime}}} = {\sum\limits_{w_{im} \in \pi_{k}}{\alpha_{nr}H_{ki}w_{im}}}}{{w_{{jm}^{\prime}} \in \pi_{k}^{\prime}},{k = 1},2,\cdots \mspace{14mu},K}} & (9) \end{matrix}$

In which, r=1, 2, . . . , N−L_(k), n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, and α_(nr) represents a coefficient of the rth interference vector in the nth vector equation. Obviously, α_(nr) cannot always be zero in each equation. Particularly, if L_(k)=N, the kth receiving end has no interference space, and M×1 zero vectors shall be put into π_(k).

III. The process of step 102 is specifically described as follows.

For the convenience of subsequent calculation of the precoding matrix, Formula (8) is converted to be:

$\begin{matrix} {{{\alpha_{n,0}H_{kj}w_{{jm}^{\prime}}} + {\sum\limits_{w_{im} \in \pi_{k}}{\alpha_{nr}H_{ki}w_{im}}}} = 0} & (11) \end{matrix}$

Formula (II) is tenable to each W_(jm′)επ′_(k), wherein k=1, 2, . . . , K, r=1, 2, . . . , N−L_(k), n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, and α_(n,0)=−1, thus Formula (II) can be represented in the form of matrix:

$\begin{matrix} {\begin{bmatrix} {{\overset{\_}{H}}_{1}\alpha_{1}} \\ {{\overset{\_}{H}}_{2}\alpha_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - N})}K}\alpha_{{({L - N})}K}} \end{bmatrix} = 0} & (12) \end{matrix}$

In which, α_(n)=└α_(n,0) α_(n,2) . . . α_(n,(N-L) _(k) ₎┘; H _(n) is an N×(N−L_(k)+1) matrix in which the first column is H_(kj)w_(jm′), w_(jm′)επ′_(k), and the others are N−L_(k) interference vectors in π_(k), i.e., H_(ki)w_(im′) w_(im)επ_(k).

Formula (11) may be further rearranged as,

Aw=0

In which, w=[w₁₁ ^(T) W₁₂ ^(T) . . . W_(1,L) ₁ ^(T) . . . w_(K1) ^(T) w_(K2) ^(T) . . . w_(K,L) _(K) ^(T)]^(T) is an LM×1 cascade precoding vector, wherein L1 is the number of the signal flows received by the first UE, and L_(K) is the number of the signal flows received by the Kth UE; correspondingly, there are L precoding vectors, and A is an (L−N)NK×LM matrix:

$\begin{matrix} {A = \begin{bmatrix} A_{11} & A_{12} & \cdots & A_{1L} \\ A_{21} & A_{22} & \cdots & A_{2L} \\ \vdots & \vdots & \ddots & \vdots \\ A_{{{({L - N})}K},1} & A_{{{({L - N})}K},2} & \cdots & A_{{{({L - N})}K},L} \end{bmatrix}} & (13) \end{matrix}$

In which, each submatrix A_(nl) is N×M dimensional, n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, l=1, 2, . . . , L and k=1, 2, . . . , K. Since each w corresponds to one data flow l, the A_(nl) corresponding to a remaining interference vector in the nth equation is a product of α_(n,0) and a channel matrix in the remaining interference vector, while the A_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector. Specifically, corresponding A_(nl) in the nth equation may be represented as

$\begin{matrix} {A_{nl} = \left\{ \begin{matrix} {\alpha_{n\; 0}H_{kj}} & {\left. l\rightarrow w_{{jm}^{\prime}} \right. \in \pi_{k}^{\prime}} \\ {\alpha_{n\; r}H_{ki}} & {\left. l\rightarrow w_{im} \right. \in \pi_{k}} \\ 0 & {others} \end{matrix}\quad \right.} & (14) \end{matrix}$

In which, l→w_(jm′)επ′_(k) indicates that in the nth equation, a vector corresponding to the lth data flow is a remaining interference vector, i.e., the vector corresponding to the lth data flow is W_(jm′)επ′_(k); l→w_(im)επ_(k) indicates that in the nth equation, a vector corresponding to the lth data flow is an interference space basis vector, i.e., the vector corresponding to the lth data flow is W_(im)επ_(k).

In order to clearly describe the deformation of Formula (11), an example is given as follows: assuming that K=4, L=6, and the flow distribution is [2 2 1 1], i.e., the first and second UEs each has two flows, the third and fourth UEs each has one flow, and the number of the transmitting antennas M and the number of receiving antennas N are 4, respectively. So L−N=2, and the IA equation is

$\left\{ \begin{matrix} {{{a_{10}H_{12}w_{21}} + {\alpha_{11}H_{13}w_{3}} + {\alpha_{12}H_{14}w_{4}}} = 0} \\ {{{a_{20}H_{12}w_{22}} + {\alpha_{21}H_{13}w_{3}} + {\alpha_{22}H_{14}w_{4}}} = 0} \\ {{{a_{30}H_{21}w_{11}} + {\alpha_{31}H_{23}w_{3}} + {\alpha_{32}H_{24}w_{4}}} = 0} \\ {{{a_{40}H_{21}w_{12}} + {\alpha_{41}H_{23}w_{3}} + {\alpha_{42}H_{24}w_{4}}} = 0} \\ {{{a_{50}H_{31}w_{11}} + {\alpha_{51}H_{32}w_{21}} + {\alpha_{52}H_{32}w_{22}} + {\alpha_{53}H_{34}w_{4}}} = 0} \\ {{{a_{60}H_{31}w_{12}} + {\alpha_{61}H_{32}w_{21}} + {\alpha_{62}H_{32}w_{22}} + {\alpha_{63}H_{34}w_{4}}} = 0} \\ {{{a_{70}H_{41}w_{11}} + {\alpha_{71}H_{42}w_{21}} + {\alpha_{72}H_{42}w_{22}} + {\alpha_{73}H_{43}w_{3}}} = 0} \\ {{{a_{80}H_{41}w_{12}} + {\alpha_{81}H_{42}w_{21}} + {\alpha_{82}H_{42}w_{22}} + {\alpha_{83}H_{43}w_{3}}} = 0} \end{matrix}\quad \right.$

Thus,

${{\overset{\_}{H}}_{1} = \begin{bmatrix} {H_{12}w_{21}} & {H_{13}w_{3}} & {H_{14}w_{4}} \end{bmatrix}},\cdots$ ${{\overset{\_}{H}}_{4} = \begin{bmatrix} {H_{21}w_{12}} & {H_{23}w_{3}} & {H_{24}w_{4}} \end{bmatrix}},\cdots$ ${{\overset{\_}{H}}_{5} = \begin{bmatrix} {H_{31}w_{11}} & {H_{32}w_{21}} & {H_{32}w_{22}H_{34}w_{4}} \end{bmatrix}},\cdots$ ${{\overset{\_}{H}}_{8} = \begin{bmatrix} {H_{41}w_{12}} & {H_{42}w_{21}} & {H_{42}w_{22}H_{43}w_{3}} \end{bmatrix}},{\alpha_{1} = \begin{bmatrix} \alpha_{10} \\ \alpha_{11} \\ \alpha_{12} \end{bmatrix}},\cdots \mspace{14mu},{\alpha_{4} = \begin{bmatrix} \alpha_{40} \\ \alpha_{41} \\ \alpha_{42} \end{bmatrix}},{\alpha_{5} = \begin{bmatrix} \alpha_{50} \\ \alpha_{51} \\ \alpha_{52} \\ \alpha_{53} \end{bmatrix}},\cdots \mspace{14mu},{\alpha_{8} = \begin{bmatrix} \alpha_{80} \\ \alpha_{81} \\ \alpha_{82} \\ \alpha_{83} \end{bmatrix}}$ w = [w₁₁^(T), w₁₂^(T), w₂₁^(T), w₂₂^(T), w₃^(T), w₄^(T)]^(T) $A = \begin{bmatrix} 0 & 0 & {\alpha_{10}H_{12}} & 0 & {\alpha_{11}H_{13}} & {\alpha_{12}H_{14}} \\ 0 & 0 & 0 & {\alpha_{20}H_{12}} & {\alpha_{21}H_{13}} & {\alpha_{22}H_{14}} \\ {\alpha_{30}H_{21}} & 0 & 0 & 0 & {\alpha_{31}H_{23}} & {\alpha_{32}H_{24}} \\ 0 & {\alpha_{40}H_{21}} & 0 & 0 & {\alpha_{41}H_{23}} & {\alpha_{42}H_{24}} \\ {\alpha_{50}H_{31}} & 0 & {\alpha_{51}H_{32}} & {\alpha_{52}H_{32}} & 0 & {\alpha_{53}H_{34}} \\ 0 & {\alpha_{60}H_{31}} & {\alpha_{61}H_{32}} & {\alpha_{62}H_{32}} & 0 & {\alpha_{63}H_{34}} \\ {\alpha_{70}H_{41}} & 0 & {\alpha_{71}H_{42}} & {\alpha_{72}H_{42}} & {\alpha_{73}H_{43}} & 0 \\ 0 & {\alpha_{80}H_{41}} & {\alpha_{81}H_{42}} & {\alpha_{82}H_{42}} & {\alpha_{83}H_{43}} & 0 \end{bmatrix}$

Next, the IA scheme is converted into an optimization issue:

$\begin{matrix} {J = {{\sum\limits_{n = 1}^{{({L - N})}K}{\alpha_{n}^{H}{\overset{\_}{H}}_{n}^{H}{\overset{\_}{H}}_{n}\alpha_{n}}} = {w^{H}A^{H}{Aw}}}} & (15) \end{matrix}$

In which, the constraint conditions are:

w _(ki) ^(H) w _(ki) =P ₁ , l=1,2, . . . ,L

α_(n) ^(T) b=−1, n=1,2, . . . , (L−N)K

Herein b=[1, 0, . . . , 0]^(T), and w_(ki) ^(H)W=P₁ indicates that a product of the transposition of a precoding vector corresponding to the ith data flow of the kth UE and the precoding vector is the power P₁ of the data flow, wherein k=1, 2, 3 . . . K ; i=1, 2, . . . L_(k), and P₁ is the power of the lth data flow in the whole data flow. For example, L=6 and the flow distribution is [2 2 1 1], i.e., the first and second UEs each has two flows, and the third and fourth UEs each has one flow, thus the first data flow of the second UE is the third flow in the whole data flow, and a product of the transposition of a precoding vector corresponding to the first data flow of the second UE and the precoding vector is P₃.

When J is small enough (approaching zero), it may be deemed that an optimal solution satisfying Formula (15) is found, which is also the solution for Formula (8).

Assuming that D is an L×LM block matrix:

D=diag{w ₁₁ ^(H) w ₁₂ ^(H) . . . w _(1,L) ₁ ^(H) . . . w _(K1) ^(H) w _(K2) ^(H) . . . w _(K,L) _(K) ^(H)}  (16)

Then, the power limitation may be represented as

Dw=e  (17)

In which, e=[P1, P2, . . . , PL]^(T). Thus, Formula (15) may be equivalent to a minimization of the following first Lagrange function:

$\begin{matrix} {{L_{c}\left( {w,\alpha_{n}} \right)} = {{w^{H}A^{H}{Aw}} + {\mu^{T}\left( {e - {Dw}} \right)} + {\sum\limits_{n = 1}^{{({L - N})}K}{\lambda_{n}\left( {1 + {\alpha_{n}^{H}b}} \right)}}}} & (18) \end{matrix}$

in which, λn and μ are Lagrange factors, and D is deemed as a known matrix, thus the complex differential of w is as follows:

∂L _(c) /∂w*=2A ^(H) Aw−D ^(H)μ  (19)

Assuming Formula (19)=0, then

w=(A ^(H) A)⁻¹ D ^(H) └D(A ^(H) A)⁻¹ D ^(H)┘⁻¹ e  (20)

Based on Formula (15), Formula (18) may be rearranged as the following second Lagrange function:

$\begin{matrix} {{L_{c}\left( {w,\alpha_{n}} \right)} = {{\sum\limits_{n = 1}^{{({L - N})}K}{\alpha_{n}^{H}{\overset{\_}{H}}_{n}^{H}{\overset{\_}{H}}_{n}\alpha_{n}}} + {\lambda_{n}\left( {1 + {\alpha_{n}^{H}b}} \right)} + {\mu^{H}\left( {e - {Dw}} \right)}}} & (21) \end{matrix}$

the complex differential of α_(n) is as follows:

αL _(c)/∂α*_(n)=2 H _(n) ^(H) H _(n)α_(n)+λ_(n) b  (22)

Assuming Formula (22) equal to 0, then

$\begin{matrix} {\alpha_{n} = {- \frac{\left( {{\overset{\_}{H}}_{n}^{H}{\overset{\_}{H}}_{n}} \right)^{- 1}b}{{b^{T}\left( {{\overset{\_}{H}}_{n}^{H}{\overset{\_}{H}}_{n}} \right)}^{- 1}b}}} & (23) \end{matrix}$

Formula (23) requires the matrix ( H _(n) ^(H) H _(n)) to have a full rank, otherwise α_(n) may be directly set as the zero space vector of ( H _(n) ^(H) H _(n)), e.g., the minimal eigenvalue vector, and normalized so that its first element is −1.

The preset initial value of α_(n) or the preset initial value w is substituted into Formulas (23) and (20), so as to orderly update α_(n) and w through the cyclic iterations of Formulas (23) and (20), and acquire an optimal solution when the preset convergence condition is satisfied. The preset convergence condition may be that the difference between Lc values of two continuous iterations is smaller than a preset first threshold, or that the number of iterations reaches a preset second threshold, i.e., a preset number of iterations.

Particularly, in consideration of the limitation of the total transmitting power, the optimization issue may be converted to be

$\begin{matrix} {{{\sum\limits_{k = 1}^{K}{\sum\limits_{i = 1}^{L_{k}}{w_{ki}^{H}w_{ki}}}} = P_{T}}{{{\alpha_{n}^{T}b} = {- 1}},{n = 1},2,\cdots \mspace{14mu},{\left( {L - N} \right)K}}} & (24) \end{matrix}$

In which, PT represents the total power for transmitting respective data flows at the transmitting end; w_(ki) ^(H)w_(ki) represents a product of the transposition of a precoding vector corresponding to the ith data flow of the kth UE and the precoding vector; and

${\sum\limits_{k = 1}^{K}{\sum\limits_{i = 1}^{L_{k}}{w_{ki}^{H}w_{ki}}}} = P_{T}$

represents that a sum of the products of the transposition of a precoding vector corresponding to each data flow of each receiving end and the precoding vector is the total power of the respective data flows.

It is easily acquired that when α_(n) is given, the optimal w is:

$\begin{matrix} {w = {\frac{1}{\sqrt{P_{T}}}u_{m}}} & (25) \end{matrix}$

In which, u_(m) is a unit eigenvector corresponding to the minimum eigenvector A^(H)A. Through the cyclic iterations of Formulas (23) and (25), α_(n) and w may be orderly updated to acquire the optimal solution.

The solution of the precoding matrix through iterations is briefly described as follows. For example, the update of α_(n) and w is regarded as an iteration cycle. It is clear that in the ith iteration, for any α^((i−1)) value, w_(n) ^((i)) value may be acquired through Formula (20). Then α_(n) ^((i−1)) value and w^((i)) value are substituted into the Lagrange function (18) while w_(n) ^((i)) value is substituted into Formula (23) to acquire α^((i)) value. After that, α_(n) ^((i)) value and w^((i)) value are used for the next iteration cycle. Thus, the following is obtained:

. . . ≧L _(c)(α_(n) ^((i−1)) ,w ^((i−1)))≧L _(c)(α_(n) ^((i−1)) ,w ^((i)))≧L _(c)(α_(n) ^((i)) ,w ^((i)))≧L _(c)(α_(n) ^((i)) ,w ^((i+1)))≧L _(c)(α_(n) ^((i+1)) ,w ^((i+1)))≧ . . .

This indicates that the algorithm is monotone but not increasing, and tends to be converged to a small value. Thus the above convergence condition can be employed.

The aforementioned IA may be performed in an interference space of either the receiving end or the transmitting end.

Specifically, for the kth UE, there are L−L_(k) interference vectors in total:

H _(ik) ^(H) g _(im) , m=1,2, . . . , L _(i) , i=1,2, . . . ,K, i≠k

when it is required to align the L−L_(k) interference vectors into an interference space of M−L_(k) dimensions,

{H _(jk) ^(H) g _(jm′) :g _(jm′)εφ′_(k)}⊂span[H _(ik) ^(H) g _(im) :g _(im)εφ_(k) ], k=1,2, . . . ,K

Herein, the definitions of φ_(k) and φ′_(k) are the similar to those of π_(k) and π′_(k), and herein are omitted. Similarly, the following IA equation may be acquired:

${H_{jk}^{H}g_{{jm}^{\prime}}} = {\sum\limits_{g_{im} \in \phi_{k}}{\beta_{nr}H_{ik}^{H}g_{im}}}$ g_(jm^(′)) ∈ ϕ_(k)^(′), k = 1, 2, ⋯  , K

In which, r=1, 2, . . . , M−L_(k), n=(k−1)(L−M)+n′, n′=1,2, . . . , L−M, and β_(nr) represents a coefficient of the rth interference vector in the nth vector equation. Obviously, βnr cannot always be zero in each equation. Similarly, if L_(k)=M, there is no interference space. In that case, N×1 zero vectors may be taken as the elements of φ_(k). Thus the IA method is still applicable.

When a set of all the equilibrium vectors of all the UEs is obtained, i.e., g_(k)=[g_(k1) g_(k2) . . . g_(kL) _(k) ]^(T) (k=1, 2, . . . , K), a precoding vector w_(k) of each of the UEs may be acquired through the following equation:

$w_{k} = {\left( {{\sum\limits_{i = 1}^{K}{H_{ik}^{H}g_{i}g_{i}^{H}H_{ik}}} + {{ɛ\sigma}^{2}I_{M}}} \right)^{- 1}H_{kk}^{H}g_{k}}$

In which, ε=0 indicates that a zero-forcing equilibrium method is employed, ε=1 indicates that a minimum mean square error equilibrium method is employed, and σ² is the noise power. Of course, w_(k) may be acquired through g_(k) in other method, without affecting the implementation of the present invention.

FIG. 3 is a schematic diagram of a digital simulation result of spectrum efficiencies acquired according to the above embodiment of the present invention. In which, M is the number of the transmitting antennas, N is the number of the receiving antennas, K is the number of the UEs, L is the total number of the signal flows, and the digits in the brackets represent the numbers of the signal flows of respective UEs. It can be seen that no BER floor occurs even in case of a high SNR.

Referring to FIGS. 4A and 4B, the embodiment of the present invention provides a device for acquiring a precoding matrix, which may be a central control apparatus independent from the transmitting and receiving ends, or located at the transmitting end. Specifically, the device includes:

an acquisition unit 401 configured to select, for each receiving end, interference vectors as interference space basis vectors from whole interference vectors;

a linear combination unit 402 configured to represent remaining interference vectors as a linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors; and

a calculation unit 403 configured to determine a currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition.

In which, the interference vector may be represented as a product of the channel vector and a precoding vector.

In that case, the acquisition unit 401 is configured to determine, for the kth receiving end, that the number of the interference space basis vectors is N−L_(K), and select N−L_(K) interference vectors as the interference space basis vectors from L−L_(K) interference vectors, wherein N is the number of the receiving antennas of the kth receiving end, L_(K) is the number of the signal flows of the kth receiving end, and L is the total number of the signal flows of a plurality of receiving ends.

The calculation unit 403 includes:

a first function acquisition subunit 4031 configured to acquire a matrix

$\begin{bmatrix} {{\overset{\_}{H}}_{1}\alpha_{1}} \\ {{\overset{\_}{H}}_{2}\alpha_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - N})}K}\alpha_{{({L - N})}K}} \end{bmatrix} = 0$

according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquire a first Lagrange function in an optimization method according to the matrix and the current channel condition; wherein α_(n)=└α_(n,0) α_(n,1) α_(n,2) . . . α_(n,(N-L) _(k) ₎┘, n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, α_(n,0)=−1, α_(n,1), α_(n,2), α_(n,3) and α_(n(N-Lk)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors;

a second function acquisition subunit 4032 configured to acquire that Aw=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquire a second Lagrange function in an optimization method according to the fact that Aw=0 and the current channel condition; l≦L, L is the total number of the data flows, an A_(nl) corresponding to a remaining interference vector in the nth equation is a product of α_(n,0) and a channel matrix in the remaining interference vector, while an A_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector;

a first iteration function acquisition subunit 4033 configured to acquire a cross-correlation function between vector a_(n) and precoding matrix w using the first and second Lagrange functions;

a first iteration subunit 4034 configured to substitute a preset initial value of w or a preset initial value of vector a_(n) into the cross-correlation function to acquire a current value of w or a current value of vector a_(n), then substitute the current value of w or the current value of vector a_(n) into the cross-correlation function until a preset condition is satisfied. The preset condition is that the difference between a value acquired by substituting the acquired current value of w or the acquired current value of vector a_(n) into the first or second Lagrange function this time and a value acquired by the substitution last time reaches a first preset threshold, or that the number of iterations reaches a second threshold.

Or, the interference vector may be represented as a product of a channel vector and an equilibrium vector.

In that case, the acquisition unit 401 is configured to determine, for the kth receiving end, that the number of the interference space basis vectors is M−L_(K), and select M−L_(K) interference vectors as the interference space basis vectors from L−L_(K) interference vectors; wherein, M is the number of transmitting antennas of a transmitting end pairing with the kth receiving end, L_(K) is the number of the signal flows of the kth receiving end, and L is the total number of the signal flows of a plurality of receiving ends.

The calculation unit 403 includes:

a third function acquisition subunit 4036 configured to acquire a matrix

$\begin{bmatrix} {{\overset{\_}{H}}_{1}\beta_{1}} \\ {{\overset{\_}{H}}_{2}\beta_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - M})}K}\beta_{{({L - M})}K}} \end{bmatrix} = 0$

according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquire a third Lagrange function in an optimization method according to the matrix and the current channel condition; wherein, β_(n)=└β_(n,0) β_(n,1) β_(n,2) . . . β_(n,(M-L) _(k) ₎┘, n=(k−1)(L−M)+n′,% n′=1, 2, . . . , L−N, β_(n,0)=−1, β_(n,1), β_(n,2), β_(n,3) and β_(n(M-Lk)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors;

a fourth function acquisition subunit 4037 configured to acquire that BG=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end; and acquire a fourth Lagrange function in an optimization method according to the fact that BG=0 and the current channel condition; wherein, each submatrix B_(nl) in B is M×N dimensional, l≦L, L is the total number of the data flows; a B_(nl) corresponding to a remaining interference vector in the nth equation is a product of β_(n,0) and a channel matrix in the remaining interference vector, while a B_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector;

a second iteration function acquisition subunit 4038 configured to acquire a cross-correlation function between vector β_(n) and equilibrium matrix G using the third and fourth Lagrange functions;

a second iteration subunit 4039 configured to substitute a preset initial value of equilibrium matrix G or a preset initial value of vector β_(n) into the cross-correlation function to acquire a current value of equilibrium matrix G or a current value of vector β_(n), then substitute the current value of equilibrium matrix G or the current value of vector β_(n) into the cross-correlation function until a preset convergence condition is satisfied. The preset convergence condition is that the difference between a value acquired by substituting the acquired current value of equilibrium matrix G or the acquired current value of vector β_(n) into the third or fourth Lagrange function this time and a value acquired by the substitution last time reaches a third preset threshold, or that the number of iterations reaches a fourth threshold.

In the embodiments of the present invention, for each receiving end, the interference vectors to be aligned are all represented as a linear combination of the interference space basis vectors, thus the interferences can be completely aligned and eliminated, without any BER floor occurring.

To be noted, in the embodiments of the present invention, each pair of transmitting end and corresponding receiving end may have different antenna configurations, without affecting the implementation of the present invention.

A person skilled in the art shall be appreciated that all or a part of steps for implementing the above method embodiments may be completed by instructing relevant hardware through a program that may be stored in a computer readable storage medium, such as ROM, magnetic disk, optical disk, etc.

The method for acquiring a precoding matrix and a device therefor provided by the embodiments of the present invention are detailedly introduced as above, and specific examples are given to illustrate the principle and implementations of the present invention. However, the descriptions of the embodiments are just used to promote the understanding of the method and the core idea of the present invention. Meanwhile, a person skilled in the art may change the specific implementations and the application ranges based on the ideas of the present invention. In conclusion, the content of the Specification shall not be construed as a limitation to the present invention. 

1. A method for acquiring a precoding matrix comprising: for each receiving end in a Multi-Input Multi-Output (MIMO) system, selecting interference vectors as interference space basis vectors from whole interference vectors, and representing remaining interference vectors as a linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors; and determining a currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and a current channel condition.
 2. The method according to claim 1, wherein, the interference vector is represented as a product of a channel vector and a precoding vector; for each receiving end, selecting the interference vectors as the interference space basis vectors from the whole interference vectors comprises: for the kth receiving end, determining that the number of the interference space basis vectors is N−L_(K); and selecting N−L_(K) interference vectors as the interference space basis vectors from L−L_(K) interference vectors; wherein, N is the number of receiving antennas of the kth receiving end, L_(K) is the number of signal flows of the kth receiving end, and L is the total number of signal flows of a plurality of receiving ends.
 3. The method according to claim 2, wherein, determining the currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition comprises: acquiring a matrix $\begin{bmatrix} {{\overset{\_}{H}}_{1}\alpha_{1}} \\ {{\overset{\_}{H}}_{2}\alpha_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - N})}K}\alpha_{{({L - N})}K}} \end{bmatrix} = 0$ according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquiring a first Lagrange function in an optimization method according to the matrix and the current channel condition; wherein, α_(n)=└α_(n,0) α_(n,1) α_(n,2) . . . α_(n,(N-L) _(k) ⁾┘, n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, α_(n,0)=−1, α_(n,1), α_(n,2), α_(n,3) and α_(n (N−Lk)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors; acquiring that Aw=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end; acquiring a second Lagrange function in an optimization method according to the fact that Aw=0 and the current channel condition; wherein, each submatrix A_(nl) in A is N×M dimensional, l≦L, L is the total number of data flows; an A_(nl) corresponding to a remaining interference vector in the nth equation is a product of α_(n,0) and a channel matrix in the remaining interference vector, while an A_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector; acquiring a cross-correlation function between vector a_(n) and precoding matrix w using the first and second Lagrange functions; and substituting a preset initial value of w or a preset initial value of vector a_(n) into the cross-correlation function to acquire a current value of w or a current value of vector a_(n), then substituting the current value of w or the current value of vector a_(n) into the cross-correlation function until a preset condition is satisfied; the preset condition is that the difference between a value acquired by substituting the acquired current value of w or the acquired current value of vector a_(n) into the first or second Lagrange function this time and a value acquired by the substitution last time reaches a first preset threshold, or that the number of iterations reaches a second threshold.
 4. The method according to claim 1, wherein, the interference vector is represented as a product of a channel vector and an equilibrium vector; for each receiving end, selecting the interference vectors as the interference space basis vectors from the whole interference vectors comprises: for the kth receiving end, determining that the number of the interference space basis vectors is M−L_(K); and selecting M−L_(K) interference vectors as the interference space basis vectors from L−L_(K) interference vectors; wherein, M is the number of transmitting antennas of a transmitting end pairing with the kth receiving end, L_(K) is the number of signal flows of the kth receiving end, and L is the total number of signal flows of a plurality of receiving ends.
 5. The method according to claim 4, wherein, determining the currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and the current channel condition comprises: acquiring a matrix $\begin{bmatrix} {{\overset{\_}{H}}_{1}\beta_{1}} \\ {{\overset{\_}{H}}_{2}\beta_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - M})}K}\beta_{{({L - M})}K}} \end{bmatrix} = 0$ according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquiring a third Lagrange function in an optimization method according to the matrix and the current channel condition; wherein, β_(n)=└β_(n,0) β_(n,1) β_(n,2) . . . β_(n,(M-L) _(k) ₎┘, n=(k−1)(L−M)+n′, n′=1, 2, . . . , L−N, β_(n,0)=−1, β_(n,1), β_(n,2), β_(n,3) and β_(n(M-Lk)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors; acquiring that BG=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end; and acquiring a fourth Lagrange function in an optimization method according to the fact that BG=0 and the current channel condition; wherein, each submatrix B_(nl) in B is M×N dimensional, l≦L, L is the total number of the data flows; a B_(nl) corresponding to a remaining interference vector in the nth equation is a product of β_(n,0) and a channel matrix in the remaining interference vector, while a B_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector; acquiring a cross-correlation function between vector β_(n) and equilibrium matrix G using the third and fourth Lagrange functions; substituting a preset initial value of equilibrium matrix G or a preset initial value of vector β_(n) into the cross-correlation function to acquire a current value of equilibrium matrix G or a current value of vector β_(n), then substituting the current value of equilibrium matrix G or the current value of vector β_(n) into the cross-correlation function until a preset convergence condition is satisfied; the preset convergence condition is that the difference between a value acquired by substituting the acquired current value of equilibrium matrix G or the acquired current value of vector β_(n) into the third or fourth Lagrange function this time and a value acquired by the substitution last time reaches a third preset threshold, or that the number of iterations reaches a fourth threshold.
 6. The method according to claim 3, wherein, a constraint condition of the optimization method comprises that the product of the transposition of a precoding vector corresponding to each data flow of each receiving end and the precoding vector is the power of the data flow; or a constraint condition of the optimization method comprises that a sum of the products of the transposition of a precoding vector corresponding to each data flow of each receiving end and the precoding vector is a total power of respective data flows.
 7. A device for acquiring a precoding matrix comprising: an acquisition unit configured to select, for each receiving end, interference vectors as interference space basis vectors from whole interference vectors; a linear combination unit configured to represent remaining interference vectors as a linear combination of the interference space basis vectors, wherein the remaining interference vectors are the interference vectors in the whole interference vectors of the receiving end except the interference vectors which act as the interference space basis vectors; and a calculation unit configured to determine a currently used precoding matrix according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end and a current channel condition.
 8. The device according to claim 7, wherein, the interference vector is represented as a product of a channel vector and a precoding vector; the acquisition unit is configured to determine, for the kth receiving end, that the number of the interference space basis vectors is N−L_(K), and select N−L_(K) interference vectors as the interference space basis vectors from L−L_(K) interference vectors, wherein N is the number of the receiving antennas of the kth receiving end, L_(K) is the number of signal flows of the kth receiving end, and L is the total number of signal flows of a plurality of receiving ends.
 9. The device according to claim 8, wherein, the calculation unit comprises: a first function acquisition subunit configured to acquire a matrix $\begin{bmatrix} {{\overset{\_}{H}}_{1}\alpha_{1}} \\ {{\overset{\_}{H}}_{2}\alpha_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - N})}K}\alpha_{{({L - N})}K}} \end{bmatrix} = 0$ according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquire a first Lagrange function in an optimization method according to the matrix and the current channel condition; wherein α_(n)=└α_(n,0) α_(n,1) α_(n,2) . . . α_(n,(N-L) _(k) ₎┘, n=(k−1)(L−N)+n′, n′=1, 2, . . . , L−N, α_(n,0)=−1, α_(n,1), α_(n,2), α_(n,3) and α_(n(N-Lk)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors; a second function acquisition subunit configured to acquire that Aw=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquire a second Lagrange function in an optimization method according to the fact that Aw=0 and the current channel condition; wherein each submatrix A_(nl) in A is N×M dimensional, l≦L, L is the total number of the data flows; an A_(nl) corresponding to a remaining interference vector in the nth equation is a product of α_(n,0) and a channel matrix in the remaining interference vector, while an A_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector; a first iteration function acquisition subunit configured to acquire a cross-correlation function between vector a_(n) and precoding matrix w using the first and second Lagrange functions; and a first iteration subunit configured to substitute a preset initial value of w or a preset initial value of vector a_(n) into the cross-correlation function to acquire a current value of w or a current value of vector a_(n), then substitute the current value of w or the current value of vector a_(n) into the cross-correlation function until a preset condition is satisfied; the preset condition is that the difference between a value acquired by substituting the acquired current value of w or the acquired current value of vector a_(n) into the first or second Lagrange function this time and a value acquired by the substitution last time reaches a first preset threshold, or that the number of iterations reaches a second threshold.
 10. The device according to claim 7, wherein, the interference vector is represented as a product of a channel vector and an equilibrium vector; the acquisition unit is configured to determine, for the kth receiving end, that the number of the interference space basis vectors is M−L_(K), and select M−L_(K) interference vectors as the interference space basis vectors from L−L_(K) interference vectors; wherein, M is the number of transmitting antennas of a transmitting end pairing with the kth receiving end, L_(K) is the number of signal flows of the kth receiving end, and L is the total number of signal flows of a plurality of receiving ends.
 11. The device according to claim 10, wherein, the calculation unit comprises: a third function acquisition subunit configured to acquire a matrix $\begin{bmatrix} {{\overset{\_}{H}}_{1}\beta_{1}} \\ {{\overset{\_}{H}}_{2}\beta_{2}} \\ \vdots \\ {{\overset{\_}{H}}_{{({L - M})}K}\beta_{{({L - M})}K}} \end{bmatrix} = 0$ according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end, and acquire a third Lagrange function in an optimization method according to the matrix and the current channel condition; wherein, β_(n)=└β_(n,0) β_(n,1) β_(n,2) . . . β_(n,(M-L) _(k) ₎┘, n=(k−1)(L−M)+n′, n′=1, 2, . . . , L−N, β_(n,0)=−1, β_(n,1), β_(n,2), β_(n,3) and β_(n(M-Lk)) are coefficients of the interference space basis vectors, respectively, when the n′th remaining interference vector of the kth receiving end in the nth vector equation represents the linear combination of the interference space basis vectors; and in H _(n), the first column is a remaining interference vector, while other columns are the interference space basis vectors; a fourth function acquisition subunit configured to acquire that BG=0 according to the linear combination of the interference space basis vectors represented by the remaining interference vectors of each receiving end; and acquire a fourth Lagrange function in an optimization method according to the fact that BG=0 and the current channel condition; wherein, each submatrix B_(nl) in B is M×N dimensional, l≦L, L is the total number of the data flows; a B_(nl) corresponding to a remaining interference vector in the nth equation is a product of β_(n,0) and a channel matrix in the remaining interference vector, while a B_(nl) corresponding to an interference space basis vector in the nth equation is a product of a channel matrix in the interference space basis vector and a coefficient of the interference space basis vector; a second iteration function acquisition subunit configured to acquire a cross-correlation function between vector β_(n) and equilibrium matrix G using the third and fourth Lagrange functions; and a second iteration subunit configured to substitute a preset initial value of equilibrium matrix G or a preset initial value of vector β_(n) into the cross-correlation function to acquire a current value of equilibrium matrix G or a current value of vector β_(n), then substitute the current value of equilibrium matrix G or the current value of vector β_(n) into the cross-correlation function until a preset convergence condition is satisfied; the preset convergence condition is that the difference between a value acquired by substituting the acquired current value of equilibrium matrix G or the acquired current value of vector β_(n) into the third or fourth Lagrange function this time and a value acquired by the substitution last time reaches a third preset threshold, or that the number of iterations reaches a fourth threshold. 